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Abstract. Neutrino and antineutrino fluxes from a core-collapse galactic supernova 
are studied, within a representative three-flavor scenario with inverted mass hierarchy 
and tiny 1-3 mixing. The initial flavor evolution is dominated by collective self- 
interaction effects, which are computed in a full three-family framework along an 
averaged radial trajectory. During the whole time span considered (t — 1-20 s), 
neutrino and antineutrino spectral splits emerge as dominant features in the energy 
domain for the final, observable fluxes. The main results can be useful for SN event 
rate simulations in specific detectors. Some minor or unobscrvablc three-family features 
(c.g, related to the muonic-tauonic flavor sector), as well as observable effects due to 
variations in the spectral input, arc also discussed for completeness. 
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1. Introduction 

The fluxes of v and V from a core-collapse supernova (SN) can encode very interesting 
information about both the (anti)neutrino properties [1] and the SN explosion 
mechanism [2]. Assuming that the latter is understood, one may try to detect distinctive 
features of flavor change in either the energy (E) or time (t) spectra of the observable 
fluxes F a (E, t), as compared with the "featureless," unoscillated fluxes F®(E, t), for 
one or more neutrino species u a . 

After the recent, seminal work in [31 H], it has been fully realized that flavor 
change phenomena driven by (anti) neutrino self-interactions [5] may induce dramatic, 
observable effects in a variety of SN scenarios, especially for inverted mass hierarchy 
[6]; see [TJ [7] for reviews of this rapidly growing research field. In particular, the so- 
called spectral split [8] or swap [9] emerges as a distinctive feature in the neutrino pjj] 
and possibly the antineutrino [TOj [11] energy spectra. If this feature dominates, then a 
stepwise flavor conversion of v e develops across a certain critical energy E c [12] , namely, 

p (E)-l 1 ( for ^<^)' m 

_ | Q (for E > ^ ^ [L) 

in terms of the survival probability P ee = P[y e - > v e) a t the end of collective effects. A 
somewhat similar feature has been observed in the V e sector [TOT ITT] , although with at 
different (generally smaller) split energy E c , 

P (E)~[ l ( E< ^)' (2) 

In a strictly 2v framework with flavors "e" and "x", the net result of collective effects 
would be a complete interchange, or "swap," of the fluxes F a between (v e , u x ) above E c 
[U, [9], [10] , as well as between (V e , V x ) above E c [TTj . 

Analogously, in a three-flavor scenario with both and u T acting as a flavor "x," 

2F X =F^ + F T , (3) 

2F X = F^ + F T , (4) 

and with usual initial conditions 

^ = 2^ = 2^ = ^=^ = ^, (5) 
a swap is expected to occur between the e flavor and one of the two x flavors, namely, 

pi = pop + pO(l _ P ) ~ J F e ( E < E cj 

(E > E N 



e 




2 K - - P«) + + F^{ 2 F f + F „ f E < g ■ (8) 

27^ =[^(l-P«) + F° I P«] + 7* I ^<^' (9) 
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In the above equations, the primes denote fluxes at the end of collective effects within 
the SN, to be further evolved up to the exit from the SN and to the final detector. 

It is important to test the above 3v expectations within comprehensive three-flavor 
calculations, for various reasons. First, one has to show that the assumed dominance of 
spectral split phenomena can indeed be manifest (without being disrupted by ordinary 
matter effects) in a sufficiently general, uncontrived and interesting supernova 3v 
scenario. Second, the effective factorization of one family out of three self-interacting 
families (leading to nonlinear equations) is never totally obvious, and "effective 2z/" 
expectations like those in Eqs. need to be explicitly checked. Third, in the self- 

interaction context, the probabilities P a p depend, among other things, on the initial 
conditions and on the absolute fluxes — which do change if the SN energy luminosity 
(which also varies in time) is distributed over 3v rather than 2v. Finally, while the v 
spectral split phenomenon is robust and largely understood in terms of lepton number 
conservation and adiabatic flavor evolution [T3J [H] , the V split seems to be more fragile 
and related to (not completely understood) nonadiabatic aspects of the evolution [TT] ; 
it is thus worth checking its appearance in a full 3v calculation. 

Several recent works have focussed on collective effects in 3v scenarios [H3, dHJ [T7J 
[T8l [T9l [20| |2T| [22] and have successfully recovered spectral split features in the neutrino 
sector (and perhaps also in the antineutrino sector [T7J [22]) in inverted hierarchy. We 
think it useful to add an independent contribution to this very recent research field, 
by discussing a SN 3u scenario where both the v and V split features are shown to 
emerge clearly for a relatively long time of ~ 20 seconds after SN explosion. Such a 
scenario is thus particularly suited to prospective experimental tests in high-statistics, 
time-integrated energy spectra of events from the next galactic core-collapse SN. 

Our work is structured as follows. In Sec. 2 we describe a representative SN 
neutrino framework, where collective effects (in the form of spectral splits) are expected 
to dominate over ordinary matter effects. In Sec. 3 we discuss the 3v formalism used 
to compute the flavor evolution of the initial fluxes F®. In Sec. 4 we present our results 
for the fluxes F' a at the end of collective effects, and show that they confirm the simple 
expectations in Eqs. ©-Q; we also discuss minor or unobservable features of our 
calculations. In Sec. 5 we complete the flavor evolution of the observable fluxes F a , by 
accounting for final phase- averaging and matter effects. In Sec. 6 we discuss the effects 
of some variations in the input SN neutrino spectra, with respect to our default scenario. 
Section 7 summarizes our results. 

The reader is reminded that our results, as well as many other observable features 
discussed in the growing literature of self-interaction neutrino effects in supernovae, 
must be taken with a grain of salt. Indeed, many open questions are still open 
from the viewpoint of the theory (validity of the mean-field approximation in the 
neutrino evolution equations), of the SN explosion energetics (neutrino luminosities and 
energy spectra) and geometry (asymmetries, turbulence), and of the neutrino evolution 
numerics (convergence and robustness of calculations). Results which currently appear 
to be rather generic (such as the spectral split phenomena in inverse hierarchy) might 
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be unpredictably challenged by improvements and further understanding in any of the 
above issues. Hopefully, the next galactic supernova explosion(s) will help to reduce 
the current level of uncertainty in the physics and astrophysics hidden in the expected 
neutrino signal. 

2. Reference neutrino parameters and SN scenario 

The adopted 3z/ oscillation parameters and our reference SN model are discussed below. 



2.1. Neutrino parameters 

We assume that the hierarchy of neutrino masses is inverted (m 3 < m 12 ), 

« ml ml) = n*±M + -Arf) , (10) 

with squared mass splittings set at the representative values 

Am 2 = 2 x 10~ 3 eV 2 , (11) 
5m 2 = 8 x 1(T 5 eV 2 . (12) 

The associated "high" (H) and "low" (L) vacuum oscillation frequencies are then 

Am 2 5.07 -,. 

u H = = — km" 1 , 13 

H 2E E/MeV [ J ' V ; 

5m 2 0.20 n u 

u L = = — ; km" 1 . (14) 

L 2E E/MeV 1 J V ; 

Within the usual parametrization [23] for the mixing matrix U = 
U(6 l2 ,6 13 , 6 23 ,5 CP ),weZx6 12 as M 

sin 2 6i2 = 0.314 , (15) 

while we consider three representative values for the more uncertain angle 823, 
corresponding to maximal and nonmaximal (but octant-symmetric) mixing: 

sin 2 9 23 = 0.50, 0.36, 0.64 . (16) 

Concerning the third angle 813, we assume a tiny reference value, 

sin 2 # 13 = 10- 6 , (17) 

in order to suppress the impact of ordinary matter effects, as explained in the following 
Subsection 2.2. Finally, we ignore possible CP-violating effects (which are arguably very 
small [22]) by setting 

S C p = . (18) 
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Figure 1. Unoscillated fluxes of neutrinos (F®, top) and antineutrinos (F a , bottom) 
with electron flavor (left) and x flavor (right). In the latter case, we plot the sum of 
muonic and tauonic fluxes. All fluxes refer to our reference supernova model and to a 
distance d = 10 kpc. The black, green, blue and magenta curves correspond to t = 1, 
5, 10, and 20 s, respectively. 



2.2. Supernova model 

Our reference SN scenario is essentially taken from [25], with minor changes in the 
average neutrino energies. We assume a galactic core-collapse supernova releasing a 
binding energy E B = 3 x 10 53 erg, equally distributed among the six (3v + 3F) species, 
and distributed in time with a time constant r = 3 s. The energy luminosity associated 
to each species is thus 

L{t) = . (19) 

6 T 

At the conventional distance for a galactic supernova, 

d = 10 kpc , (20) 
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the unoscillated flux of the neutrino species v w per unit of area, time, and energy, is 
expected to be 

F o (F f] _m €(e) , , 

where we assume normalized thermal energy spectra 4>° a {E) with average energy (E a ). 
As in [10], inspired by [26] . we set 

(E e ) = 10 MeV, (E e ) = 15 MeV, (E x ) = 24 MeV , (22) 

corresponding to temperatures (T e , T g , T x ) = (3.17, 4.76, 7.62) MevJ|| Figure 1 shows 
the corresponding fluxes F®(E, t) in the energy interval E G [1,50] MeV, at four 
representative times (t = 1, 5, 10, and 20 s). 

Concerning the SN geometry, we adopt the spherically-symmetric bulb model pQ, 
with a neutrinosphere radius R u = 10 km. Within this model, the previous initial 
conditions for the neutrino luminosity and unoscillated fluxes define, at any radius 
r > R u , the effective density of neutrinos (N = N e + + N T ) and of antineutrinos (N) 
per unit volume, as well as the self-interaction potential 

/x(r) = V2G F [N(r) + N(r)} . (23) 

For each species u a , it is useful to introduce also density n a per unit of volume and 
energy (see [4J, [TO] for details): 

N a = J dEn a (E) . (24) 

Figure 2 (left panel) shows the function fi(r), up to r = 500 km, for the four t 
values chosen. In the same panel, the shaded horizontal band corresponds to the range 
A* ^ [fMni, /4up] — [6-7, 67] km -1 where, according to the estimate in pjj], collective 
bipolar {27] ("pendulum" [6]) oscillations are expected to develop (after a "synchronized" 
regime [271 El I2H]) for this SN scenario in the 2u case. We have explicitly verified that 
these expectations also hold in our full 3u flavor evolution, with a radial accuracy better 
than ~ 10 km (not shown). For example, at t = 10 s, we find numerically a bipolar 
development range r ~ [60, 100] km, in good agreement with the range r ~ [55, 95] km 
derived from the intersection of the horizontal band with the /i(r) curve at t — 10 s. 

However, collective phenomena extend somewhat beyond the bipolar range, and 
eventually vanish when the spectral splits are complete. As a rule of thumb, the end of 
all collective effects occurs at a radius r en( j where the self-interaction potential \x equals a 
relatively small value // e nd- 111 our numerical experiments, we find that /i cn d — 0.5 km -1 
provides a reasonable estimate for the radius r em j 

/i(r end ) ~ 0.5 km" 1 , (25) 

at any time t (e.g., r en d — 400 km at t = 1 s in Fig. 2). For safety, we have always 
checked that our results do not appreciably change by continuing the numerical flavor 
evolution for further ~100 km or more (e.g., up to r ~ 500 km for t — 1 s). 



X Effects of variations with respect to the default choice in Eq. (|22|) will be discussed in Sec. 6 
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Figure 2. Radial profiles of the self-interaction potential fj, (left) and of the matter 
potential A (right). The black, green, blue and magenta curves correspond to t = 1, 
5, 10, and 20 s, respectively. In the left panel, the shaded horizontal band marks the 
H range where bipolar effects develop. In the right panel, the shaded band marks the 
range of the high vacuum frequency u>h for E € [1, 50] MeV, where MSW effects may 
develop if the iJ-resonance condition (A ~ ojh) is satisfied. 




Figure 2 (right panel) shows the radial profile of the matter potential A at different 
times, as taken from [25]. As usual, the matter potential is defined as 

X(r) = V2G F N e -(r) , (26) 

where N e -(r) is the net electron density. In inverted hierarchy, antineutrinos may 
undergo significant matter effects (or MSW effects, from Mikheev, Smirnov, and 
Wolfenstein [29|) when the so-called if-resonance condition is fulfilled (see, e.g., [30]). 

A ~ u H . (27) 

For the energy range in Fig. 1 (E e [1, 50] MeV), the relevant interval for the high 
vacuum frequency, ujh £ [0.1, 5] km -1 , is shown in the right panel of Fig. 2 as a shaded 
horizontal band. [Note in Fig. 2 that the uoh value of 0.5 km -1 in Eq. ( 125]) is in the 
middle of such band.] Its intersection with one of the A(r) curves marks the related 
if -resonance range (for example, from Fig. 2, r > 130 km at t = 10 s). 

A comparison of the left and right panels in Fig. 2 shows that, in our reference SN 
model, collective bipolar oscillations develop before possible if-resonance effects at any 
time (while they may be coupled in other scenarios with shallower matter profiles, see 
e.g. [HUH]). Flavor swap effects, however, extend up to r cnd estimated previously, and 
may partly enter the if-resonance region for large ujh at t > 1 . However, for values of 6*13 
(at least) as small as in Eq. ffTT|) . the H resonance is highly nonadiabatic, and produces 
no further flavor conversion at any relevant SN energy [30]. We have numerically verified 
the absence of MSW flavor conversion for antineutrinos in inverted hierarchy (besides 



Supernova neutrino three-flavor evolution with dominant collective effects 



8 



their collective flavor swap) for r well within the if -resonance region. We have also 
verified that all our main results in Sees. 4 and 5 remain unchanged for reasonably 
smaller values of sin 2 #13 (e.g., 10~ 7 or 10~ 8 ), which simply cause a logarithmic delay of 
the bipolar oscillation onset [21 [6] . As an added bonus, for sin 2 13 < 1CT 6 , one gets a 
basically complete suppression of: (i) possible matter effects due to shock-wave features 
(see, e.g., [31]); {ii) possible decoherence effects due to density fluctuations (see, e.g., 
[32]); (Hi) Earth matter effects, if any, before detection (see, e.g., [TJ,[20]). These welcome 
simplifications make neutrino-neutrino interaction effects dominate over neutrino-matter 
interaction effects at any time t > 1 s in our framework, and motivate a posteriori our 
choice for the tiny mixing angle in Eq. (I17jl . 

A final remark is in order. As described in the next Sections, we evolve the 
neutrino flavors along a single, averaged radial trajectory. This (so-called "single-angle") 
approximation is often reasonable — when compared with few available "multi-angle" 
calculations — but it may fail badly when the matter potential is so high that its gradient 
depends sensitively on the trajectory (see [33J, and references therein). For our inspiring 
SN model [25], this situation may occur only at very early stages (t < 1 s), which are 
excluded from the present investigation. 



3. Three-flavor formalism 



The density matrix formalism is particularly useful to deal with neutrino self- 
interactions, which depend on the neutrino density itself. In the two-family case, the 
traceless part of the (2 x 2) density matrix p can be expanded onto Pauli matrices <7j, 
with projections coefficients forming a 3- vector P for each energy mode E [6]. 

In the three-family case, we have generalized the 1v evolution equations to the 2>v 
case [31], by projecting the (3 x 3) density matrix p in flavor basis onto Gell-Mann 
matrices Aj, in terms of an 8-component vector P. The resulting equations [M] are 
formally similar to those recently discussed in [T71 [T8l [19], and are briefly reviewed 
below for completeness. 

Given an 8-dimensional orthonormal basis (ei,...,eg) for the P vector space, 
the flavor projector matrices, II e = diag(l, 0, 0), II M = diag(0, 1, 0) and II r = 
diag(0, 0, 1), can be represented as 

n Q = ^ + u a • ^ (a = e, ft, t) , (28) 

where 1 is the 3x3 unit matrix, A = J2 A« e^, and the Aj's are the Gell-Mann matrices 

with usual conventions: Aj = Aj; tr(Aj) = 0; tr(AjAj) = 25^; and [A*, Xj] = 2ifijk\k- 

The u Q vectors (with |u a | = 2/v^3) read: 

1 12 
u e = e 3 + -^=e 8 , u M = -e 3 + -^=e 8 , u T = -^ e 8 • (29) 

The neutrino density matrix, at given energy E, is projected via 

p = n (± + p.±) , (30) 
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where P = E^e* (with |P| = 2/v^), and n is the neutrino density per unit volume 
and energy introduced before. At any time t and energy E, the v a content is given by 
tr(p(t)n Q ). 

The evolution of p is governed by the Liouville equation i'p = [H, p], where the 
Hamiltonian H contains vacuum, matter, and self-interaction terms. In the vacuum 
term, we include both low and high vacuum frequencies {u L and Uh, respectively), and 
assume inverted hierarchy. In the matter term, we include the u T — u fl potential difference 
at one loop [35], whose size is 5\/\ ~ 5 x 10~ 5 (see also [HI [IE]). By expanding p and H 
onto Gell-Mann matrices, evolution equations are obtained for the v flavor polarization 
vector P and, analogously, for the V vector P: 



P 
P 



f (u L B L - uj h B h ) + Av + /xD] x P , (31) 
-(u L B L - u H B H ) + Av + fjD] x P , (32) 

where the 8-component vector product is defined as (a x b); = a,j bk- 

In the above equations, the first (vacuum) terms embed the squared mass splittings 
via ul,h, and the mixing angles via effective "magnetic fields" B L H with components 



B, 



Cl3(Sl2C 2 3 + ^12^13523) \ 



|Cl 2 (3 + 3Ci3 + 3C*23 — C13C23) + |s 13 5'i2>S'23 
2Cl2'S'l3 c 23 — 5'l2Cl3S23 



— S13C23S12 — — Ci3)5 , 2 3Ci2 



Ci2(l + C13 — 3C*23 + C13C23) — -^S 13 5 , i2iS'23 / 



(33) 



and 



B 



H 



i(l-3C 





3 + C23 




\ 



C13C' 



13*^23 J 



^(1 + C13) 



where we have set s 



1.1 



sin 6. 



3C 



13 





• 3C 2 3 



(34) 



3C13C23) J 



ij: Cij 



cos 6ij , Sij 



sin26ij, and 



cos 26; 



The second (matter interaction) term Av in Eqs. ( !3T|) and (|32|) embeds the v e — v^ T 
interaction energy difference (oriented along u e ), as well as the tiny correction due the 
v T — interaction energy difference (oriented along u T ), 



Av 




(35) 
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The third (self-interaction) term /xD couples all neutrino and antineutrino modes via 
the integral vector difference 



Note that, in general, neutrino self- interactions depend on the crossing angle between 
the neutrino trajectories [Sill]- We have implicitly assumed that such dependence can be 
averaged out along a single, radial trajectory (single-angle approximation, see [10J and 
references therein). If crossing angles were accounted for, the equations would entail 
a further angular variable (multi-angle description) [U [7J QUI EE] not included in the 
present investigation. 

We discretize in energy the coupled evolution equations ( 13T1) and ( l32i) . and solve 
them by numerical integration (see [10] for details), up to the end of collective effects 
(r cn( j<500 km). The initial conditions are fixed by the SN model described in the 
previous Section, for each representative time: t = 1, 5, 10 and 20 s. The results are 
described below in the graphically convenient range EE [1, 50] MeV (although the 
numerical evolution includes modes with E < 1 MeV). 

4. Intermediate fluxes after collective effects 

In this Section we discuss our numerical results for the intermediate fluxes F' a at the end 
of collective effects (r = r cn( j<500 km), and show that they confirm the expectations 
in Eqs. (JHJ) — (EI) - In order to allow a visual comparison with the unoscillated fluxes F® 
at d = 10 kpc (Fig. 1), the fluxes F' a are rescaled by a factor (r cnd /d) 2 . We stress that 
the intermediate fluxes F' a are unobservable, as both v and V are subject to further 
flavor evolution in the SN, and to phase-averaging effects up to the detector. The final, 
observable fluxes F a at d = 10 kpc will be estimated in the next Section. 

Figure 3 shows our computed fluxes F' a , in the same scale as Fig. 1. The comparison 
of Figs. 1 and 3 confirms that spectral splits of neutrinos and antineutrinos emerge as 
dominant features, as expected. In our adopted scenario, the neutrino split occurs 
around a critical energy E c ~ 7 MeV [10], dictated by lepton number conservation 
[T3| [T4"] . The antineutrino split occurs at a somewhat lower energy E c ~ few MeV, 
which is not predicted a priori so far [TT]. However, apart from a transition region 
(about 1-2 MeV wide) around E c and E c , the evolved fluxes in Fig. 3 coincide with 
those expected from the application of Eqs. (JHI) ([9]) to the unevolved fluxes of Fig. 1. 

This is shown in more detail in Fig. 4, with reference to the time t — 10 s (similar 
results hold at 1, 5, and 20 s). Below the critical energy, all u a species remain unchanged: 
the evolved fluxes F' a (solid blue curves) coincide with the unevolved fluxes F° (dashed 
blue curves) in the upper left panel. Conversely, above the critical energy, the electron 
flavor fully converts. As a consequence, in the lower left panels of Fig. 4, the evolved 
fluxes F' e and F' e turn into F° and F° x , respectively (dashed red curves). Antineutrinos 
show a similar behavior (right panels). Therefore, Eqs. ([B])-© hold with good accuracy, 




(36) 
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Figure 3. Fluxes of neutrinos (i 7 ^) and antineutrinos (F a ) at the end of collective 
effects, rescaled to d = 10 kpc for comparison with Fig. 1. See the text for details. 



except close to the spectrum step, whose details may also depend on the adopted single- 
angle approximation. Indeed, multi-angle calculations (not performed in this work) may 
further widen the transition region [TO], especially for antineutrinos [TT]. However, they 
are not expected to change the low- and high-energy limits of the spectra. 

The above Zv results, obtained at the end of collective effects, are not significantly 
influenced by the subdominant "solar" squared mass difference {5m 2 ) or by the — v T 
interaction energy difference (8 A) . Effects of 5m 2 are expected to be larger at relatively 
shallow matter densities [HUH], i.e., at larger t in Fig. 2; we find the largest fractional 
flux variations at t — 20 s, as obtained by setting 5m 2 = 0, to be negliglible (< few%). 
Conversely, effects of S X are expected to be larger at higher matter densities [T6 l [T8 l 122] . 
i.e., at shorter t; we find the largest variations at t = 1 s, as obtained by setting 
5X/X = 0, to be also negligible (< few%). We find that these variations are mainly 
localized around the critical split energies, and can thus be mainly attributed to small, 
subleading 3z/ collective effects. 

To a very good accuracy, our results are independent of any effect which may change 
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Figure 4. Fluxes of neutrinos and antineutrinos at the end of collective effects, for 
t = 10 s. Solid blue curves: computed fluxes F' a . Dashed blue and red curves: limiting 
behavior at low-energy and high-energy, respectively, in terms of unoscillated fluxes 
F°, according to Eqs. 



the relative and v T flux proportions within the sum 2F' X = F' e + F'^, including not 
only the small v^ T correction SX/X, but also the (quasi)maximal u^ T mixing angle #23- 
Indeed, the curves in Fig. 3 do not depend on the chosen value of sin 2 #23 in Eq. ( fT6l) , 
as we have explicitly checked. 

For instance, Fig. 5 shows both the absolute fluxes with e and x flavor (upper panels) 
and the muonic-to-tauonic flavor ratio (lower panels) at the end of collective effects, for 
t — 10 s; similar results hold at different t (not shown). The curves in the upper panels 
do not depend on the indicated value of sin 2 #23 , while those in the lower one do. More 
precisely, for neutrinos (lower left panel) the ratio F' JF' T remains equal to 1(= F®/F®) 
in three cases: (z) at maximal mixing, 623 = 7r/4, where and v T are interchangeable, 
up to minor b~XjX effects; (ii) below the critical energy E c ~ 7 MeV, where there is 
no net flavor conversion; and (Hi) at the specific "equalization" energy E cq ~ 20 MeV 
where F e ° = F®, and flavor conversions become inoperative: F' e = F' x . Similarly for 
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Figure 5. Neutrino and antineutrino fluxes F' a at the end of collective effects 
(r<500 km), for t = 10 s. Upper panels: absolute fluxes. Lower panels: muonic- 
to-tauonic flux ratio for v (left) and V (right), for sin 2 #23 = 0.36, 0.5, 0.64. 



antineutrinos (lower right panel), but with E c ~ 3-4 MeV and E cq ~ 27 MeV. 

In Fig. 5, the cases where F' / F' T 7^ 1 can be understood by considering that, in 
inverted hierarchy, flavor conversions occur mainly between the higher mass state z/3 
and the lower mass doublet z/ 12 , with z/ 3 being richer in v T or according to #23 being 
in the first or second octant (while the doublet is always richer in v e ). If F e ° > F®, 
namely E < E cq , then Am 2 -driven conversions are dominantly of the kind v e — > v T 
[y e —> Vp) for sin 2 ^23 < 1/2 (sin 2 #23 > 1/2), so that the final v T (z/ M ) flavor increases 
at the expenses of {y T ). The opposite happens if F e ° < F®, namely E > E eq . The 
same reasoning hold for antineutrinos. 

These arguments explain the main qualitative features of the numerical results in 
the lower panels of Fig. 5 which, by themselves, have just an academic interest. They 
may serve, however, as benchmarks in the exploration of extended SN scenarios with 
nonstandard initial conditions [37] or interactions [38] in the v^-v^- sector. 
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5. Final fluxes at the detector 

After collective effects have ended, and the if -resonance region is traversed, the flavor 
evolution is eventually subject to the so-called "L resonance" at A ~ u^, which can be 
assumed adiabatic (see, e.g., [30]). [As noted, additional effects due to possible shock 
wave features or density fluctuations are ineffective for #i 3 as low as in Eq. (TTTT) .] 

In the specific context of self-interacting SN neutrinos in three families, final L- 
resonance effects have also been numerically verified in jTTj , by continuing the numerical 
evolution up to a few thousand km. Here we do not repeat this numerical check of 
long-distance evolution, but take the resulting effects for granted; for completeness, we 
rephrase the related arguments of [17] in our notation as follows. 

In the absence of self interactions effects, a strongly nonadiabatic H resonance, 
plus an adiabatic L resonance, would eventually distribute the original fluxes over the 
effective mass eigenstates and T>i as [30] : 

Fi = , (37) 

Fi = F° e , (38) 

^3 = F° x , (39) 

Fi = F° e , (40) 

F 2 = F° X , (41) 

F 3 = F° X . (42) 

Spectral splits, however, alter these standard MSW expectations and swap the flavor 
contents in the effective (e, 3) sector [IT], so that 

(43) 
(44) 
(45) 
(46) 
(47) 
(48) 

where P ee and P ee , embedding collective effects, have a stepwise behavior as in Eqs. (pQ) 
and ([2]). The final, phase-averaged fluxes for the electron flavor at detection are given 
by 

Fe = £ \ U ei\ 2 Fi * cos 2 9 12 F l + sin 2 6 12 F 2 , (49) 
F e = J2\ U ei\ 2 F t ^ cos 2 6viF x + sin 2 6 l2 F 2 , (50) 
which, together with Eqs. (j4"3]) -(l4"8" ]) . reproduce the limiting behaviors discussed in [T7] . 
/ cos 2 6 12 F° + sin 2 6 12 F° (E < E c ) , 
6 " I ^ (E > E c ) , (51) 

, sin 2 ^ 12 F;. + cos 2 ^ (E<E C ) , ^ 

x n (F>E c) , (52) 



Fx 


= F° 




F 2 


= F e °P ee + F x °(l- 


Pee) 


F, 


= F°P ee + F e °(l- 


Pee) 


Fx 


= F° e P ee + F° x (l~ 
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Figure 6. Final oscillated fluxes of neutrinos (F a ) and antineutrinos (F a ) at 
d = 10 kpc. The comparison with the unoscillated fluxes in Fig. 1 shows the imprints 
of collective flavor evolution effects. 



where the low-energy V split, not considered in [T7j, is also included. 

By using Eqs. (ED - ©, one can eliminate P ee and P ee in terms of the fluxes F' a at 
the end of collective effects. The final fluxes F e and F e can thus be expressed in terms 



of the initial fluxes in Fig. 1 and of the intermediate computed fluxes in Fig. 3, 

F e = cos 2 6 12 F° + sin 2 6 l2 F' e , (53) 
F e = cos 2 9 12 F' e + sin 2 Q 12 ~f x . (54) 

Finally, conservation of the total v and V flux provides F x and F x by subtraction, 

2F X = (2F° X + F e °) - F e , (55) 
2F x = (2F x +7 fi e )-F e . (56) 



Figure 6 shows the corresponding results in graphical form. The spectral split 
features in the v e and V e spectra are somewhat reduced by the low-energy f9 12 mixing 
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from Eqs. (H9"j) and (1501) . but are still clearly observable at any time t > 1 s. The x-flavor 
split features are instead more suppressed by mixing. 

We emphasize that, in view of prospective observations of galactic SN neutrino 
bursts (at least for the electron flavor via charged currents) the persistence of similar 
stepwise features for several seconds is rather useful. Detection of such features (if 
realized in nature) requires setting the threshold at low energy (< few MeV): a very 
challenging goal, since one has to fight against backgrounds (which may be large in some 
shallow detector projects (39]) and low cross sections. However, low signal rates may 
be at least compensated by integration over time intervals of a few seconds, without 
canceling the persisting split effects. In order to facilitate feasibility studies in specific 
experimental settings, we provide our main results (i.e., the final fluxes at d — 10 kpc in 
Fig. 6, as compared to the unoscillated ones in Fig. 1) in computer-readable form upon 
request. 



6. Effects of variations in the input SN neutrino spectra 

Our references choice of neutrino energy spectra and luminosities is by no means unique. 
Different choices may be motivated by the results of some SN explosion simulations, 
see e.g. [ID]. We do not consider herein possible deviations from the hypothesis of 
luminosity equipartition among different neutrino species [IDJ, which are not supported 
in our reference time interval (t = 1-20 s) by the results of [26] , and which might lead to 
new and more complicated split features [JT] beyond the scope of this paper. However, 
we do study some interesting variations of the average energies (E a ). 

In particular, following Ref. [12] (inspired by [IDJ) one may adopt the values: 

(E e ) = 12 MeV, (E e ) = 15 MeV, (E x ) = 18 MeV (at any t) , (57) 



which are much closer to each other than our reference choice in Eq. ( 1221) . We shall refer 
to the choice in Eq. fl571) as to the "smaller A(E) scenario." 

One may also include average energies with noticeable changes with time, as 
suggested, e.g., by the results in [26], which read: 

(E e ) = 13 MeV, (E e ) = 16 MeV, (E x ) = 23 MeV (t = 1 s) ; (58) 

(E e ) = 11 MeV, (E e ) = 18 MeV, (E x ) = 25 MeV (t = 5 s) ; (59) 

(E e ) = 11 MeV, (E e ) = 20 MeV, (E x ) = 25 MeV (t = 10 s) ; (60) 

(E e ) = 11 MeV, (E e ) = 20 MeV, (E x ) = 25 MeV (t = 20 s) . (61) 

We shall refer to such assignments as to the "time-dependent (E) scenario." 

For fixed energy luminosity [taken as in Eq. ( 1T91) ]. variations in the average neutrino 
energies lead to variations in their total number and thus also in the self-interaction 
potential fi, as compared to the reference one shown in Fig. 2. Figure 7 shows the 
radial profile of fi for the case of smaller A(E) (left panel) and of time-dependent 
(E) (right panel). The shaded bands corresponds to the range ji G [yU inf , yU sup ] where, 
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Figure 7. Radial profiles of the self-interaction potential fj, at t = 1, 5, 10, and 20 s 
for two different variations of the supernova neutrino input (as compared with Fig. 2). 
Left: case with smaller differences among (constant) average energies. Right: case with 
time-dependent average energies. The shaded bands mark the fi range where bipolar 
effects develop. The band acquires a time dependence in the right panel. See the text 
for details. 



according to the discussion in [10] , collective bipolar oscillations occur; the range is time- 
dependent in the right panel. Both the /x(r) curves and the bands in Fig. 7 differ from the 
corresponding ones in Fig. 2. In particular, at t = 20 s, the intersection of the fi(r) curve 
with the shaded band in Fig. 7 leads to the smallest radial interval (Ar ~ 15 km) for 
the development of bipolar oscillations — a fact that leads to an interesting consequence, 
as shown below. 

Figure 8 shows the results of the neutrino flavor evolution (at the end of collective 
effects) in the case of smaller A{E). Compared with our reference scenario in Fig. 3, 
the spectral split features are qualitatively similar but less pronounced, as expected 
from the fact the spectral differences of different species are smaller. In particular, 
the step-like variation across the split is reduced by a factor of two or more, for both 
neutrinos (around E c ~ 7.5 MeV) and antineutrinos. From the experimental viewpoint, 
the smaller the split variations, the higher the energy resolution and the statistics needed 
to observe them. Future SN explosion simulations will shed new light on the expected 
(dis) similarities among the unoscillated spectra of u e , V e and u x , and thus on the size of 
observable oscillation effects (collective or not). 

Figure 9 is analogous to Fig. 8, but refers to the case of time-dependent (E) (to be 
compared with the reference scenario in Fig. 3). In this case, the critical split energies 
appear to be also time- dependent. For neutrinos, conservation of net lepton number 
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Figure 8. Fluxes of neutrinos (F' a ) and antineutrinos (F a ) at the end of collective 
effects, for the case with smaller A(E) (as compared with the reference scenario in 
Fig. 3). 



(assuming no V split in first approximation) predicts E c ~ 7.2, 8.9, 9.7 and 9.7 MeV for 
the chosen spectra at t — 1, 5, 10 and 20 s. This trend is consistent with the numerical 
results in the upper panels of Fig. 9. Antineutrinos (lower panels) also show a steadily 
increasing split energy, which we are unable to estimate a priori, however. The increase 
in critical energies with time may favorable from the experimental viewpoint, since the 
cross section increase compensate in part the luminosity decrease; on the other hand, 
it may prevent integration of spectra on time intervals larger than At ~ 1 s, where the 
split feature would be blurred. 

Interestingly, no split appear in Fig. 9 at t = 20 s: the spectra are basically 
unoscillated. This fact is related to the very short range (Ar ~ 16 km) expected 
for the development of bipolar oscillations at t = 20 s, as noted above. It turns out that 
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Fluxes of neutrinos (F' a ) and antineutrinos (F a ) at the end of collective 
the case with time-dependent (E) (as compared with the reference scenario 



the period of bipolar oscillations is of the same order of Ar in this case and, literally, the 
"flavor pendulum" has not enough time to make a single swing: it remains frozen in the 
upward, unstable equilibrium position, due to the extremely sudden decrease of /x(r). 
Conversely, in our reference scenario at t = 20 s (see Fig. 3) the range Ar happens to 
be slightly larger, thus enabling the flavor pendulum to make a full swing and to start 
the collective transitions, which eventually lead to the spectral split. At late time, the 
neutrino and antineutrino splits appears thus to be relatively fragile phenomena, whose 
absence or presence might provide, in principle, some information about the gradient of 
the self- interaction neutrino potential /x(r). 
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7. Summary 

Building upon recent literature on 3z/ collective effects in core-collapse supernovae 
[T5| [T6l [TTl [HI Uni [201 [22], we have performed an independent study of three-flavor 
evolution of neutrinos and antineutrinos within a rather "standard" SN model (Figs. 1 
and 2). Then, assuming inverted neutrino mass hierarchy and tiny 6\z (i-e., strongly 
nonadiabatic if-resonance) , self-interaction effects are expected to provide the dominant 
spectral features, in the form of "splits" — or "stepwise swaps" — for both v and V. 

We have explored this SN 3z/ scenario by solving the evolution equations for the 
2>v density matrix in vector form and single-angle approximation, after discretization in 
energy. The numerical results for the evolved fluxes at the end of collective effects 
(Fig. 3) confirm basic expectations at low and high energy [Eqs. in terms 

of unevolved fluxes (Fig. 4). Effects of the "solar" mass splitting and of the v ii -v T 
interaction energy difference in matter are found to be negligible; effects due to variations 
of 6*23 are large (Fig. 5) but unobservable. The final fluxes at the detector (Fig. 6) are 
obtained by standard application of adiabatic L-resonance effects. 

Both v and V spectral split features tipically emerge with similar characteristics 
in the whole time interval considered (t = 1-20 s). Observations of such features (if 
realized in nature) and of their possible variations require low-energy thresholds, where 
the signal is suppressed by low cross sections; however, time integration can partly 
overcome the suppression, without necessarily canceling the persistent spectral features. 
If the collected statistics and the energy resolution allow, spectra collected at different 
times might even reveal migrations of the split energies due to time variations of neutrino 
temperatures. The low-energy frontier in SN neutrino physics may thus be the key to 
access the physics of inverted mass hierarchy, especially if #13 is very small. 
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